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Abstract — This work develops novel algorithms for high- 
performance networking in the presence of obstacles based on 
a method for communicating via ultrasonic rays reflected at 
the obstacles. The rays are curves determined by the variable 
speed of sound and initial conditions and we develop ultrasonic 
ray models based on a system of differential equations. We 
present new parallel algorithms and software for shape and 
trajectory reconstruction of moving obstacles and show how the 
reconstructed reflection point of a ray at an obstacle is a natural 
router for messages between the ray's transmitter and receiver 
and discuss the advantages of the new architecture. We discuss 
how the new algorithms and software improve the performance 
and properties of the network architecture. 

I. Introduction 

Consider an ultrasonic wave, or signal, described by the 
wave equation 

utt - c^{x,y,z)Au = (1) 

where c{x, y, z) > is the variable speed of sound in 

n{t) = no\Ti[{t) and 



= 



for t > and i}i{t) C fio C M^, where ili{t) is a convex 
moving obstacle with a smooth boundary in a bounded domain 

We consider an environment without caustics and look for 
solutions of the form 



u{x,y,z,t) 



iuj{W(x,y,z) — t) 



(2) 



where the eikonal function W{x, y, z) = const defines a 
surface of constant phase. These solutions of the ray equation 
are called rays or ray solutions. Suppose that for all t we 
are given all integrals .f{l)dl = where 7 are broken 
rays in il{t) such that for each point P G d^li{t) there is 
at least one broken ray 7 reflecting at P. A broken ray is a 
ray reflecting at the obstacle and starting and ending at the 
observation boundary OHq. 

Let f(x) ~ -r-^ — r > in fi. Then as we know C~, 
correspond to signal travel times in a medium with speed 
of sound c{x,y,z). The shape and trajectory reconstruction 
problem is to find cJili (t) given the sets (t) where 7 G fl{t). 

Consider Silo as the observation boundary and signal trans- 
mitters and receivers with known locations along the boundary. 
Transmitters send ray signals with known initial incident and 



azimuth angles at known transmission times. Receivers receive 
reflected ray signals and record the time when the signal was 
received. The combined data leads to a set of data points 

Bk = {xuyi, Zi,Xr,yr, Zr,4>k,Qk,tk,(,k) 

where where (pk and Ok are the initial incident and azimuth 
angles of the ray from the transmitter, xi, yi, zi are the 
coordinates of the transmitter endpoint of the ray, and Xr, yr, 
Zr are the coordinates of the receiver endpoint of the ray, tk 
is the time of flight for the signal and £_k is a frequency of 
the signal. We present algorithms for the shape and trajectory 
reconstruction problem for finite travel times and variable 
speed of sound derived from the equations for the Shooting 
Method for two-point seismic ray tracing [JG]. 
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This system of equations derived from the eikonal equation 
has wide applications in seismology and is used in algorithms 
for seismic ray tracing [SK]. The system of equations 3 is for 
a Cartesian coordinate system where R{t) = {x{t), y{t), z{t)) 
is the ray position vector, (j){t) is the incident angle of the ray 
direction vector with the z axis and 9{t) is the azimuth angle 
that the projection of the ray direction vector makes with the 
positive X axis. 

We consider that the speed of sound c{x, y, z) is known 
inside fi. We will reconstruct the position of the reflection 
point P given the positions of a transmitter L and receiver 
S, the incident and azimuth angles of the transmitted ray 
with respect to a Cartesian coordinate system centered at the 
transmitter, and the time of flight of the signal. We know 
the initial position (x(0), y(0), z(0)) and therefore we know 
the velocity of the ray al t = 0. We know the angle and 
initial speed of wave propagation for the transmitted signal at 
time t = 0. Then, knowing the initial conditions, we can find 
numerically the signal path from the transmitter L through the 
reflection point P for a given travel time Tk- Let the travel 



time from the transmitter L to the reflection point P be and 
the total travel time from the transmitter L to the receiver S 
through the reflection point P be where tk > Tk- 

From symmetry, we can imagine that the signal received at 
the receiver S is transmitted from the receiver S and that it 
arrives at the reflection point P in time t = tk — Tk because 
its travel time is tk — Tfc. The ray path from S to P is again 
described by the above system of equations for seismic ray 
tracing, however we do not know the initial angles for the 
signal from the receiver S. 

In order to reconstruct the intersection point Pk of rays 
starting from the transmitter and receiver for a given data point 

Bk = {xuyi,Zi,Xr,Vr,Zr,(t)k,Ok,tk,ik) 

containing the measured values for the signal, we step through 
a discrete set of values for Tk- Starting to trace a ray path 
from L, at each time step, each new point on the path is a 
candidate reflection point P. At this stage, this is basic initial 
value ray tracing. Then we apply the two point seismic ray 
tracing shooting method to see if we can reach from receiver 
S candidate point P for the remainder of our time budget. If 
we can then we have found the reflection point P for this data 
point. Otherwise, we continue with the next time step of the 
initial value ray tracing and check the next candidate point P. 
We repeat this until we find a reflection point P or exhaust 
our time budget tk- For example, let 

_ tkns 

where Nr is an integer that specifies the time-step resolution 
and Us is an integer that specifies the number of the time step 
and such that Q < Us < Nr- Each time Tks corresponds to a 
unique point Pks on the ray path starting from the transmitter 
L such that Pks can be reached from L in time Tks - Next, for 
each time step corresponding to time Tks to reach candidate 
Pks from L, we step through the range of initial angles for 
the signal starting from the receiver S and check whether the 
curve from S will intersect Pks in time tk — Tks- Alternatively, 
we can step through a range of reflection angles at the point 
Pks that is found for the given Tks and check whether the 
reflected signal will intersect S in time tk — Tks- 

II. Reconstruction Algorithms 

The input to the following algorithm is the speed of sound 
c{x) for the domain Q, and a set of data points or ray 
coordinates corresponding to broken rays. The output is a 
unique set of points in R'^. The points from the output are 
the reflection points reconstructed from the input data. 
Require: Set of broken ray data points Bk = 

{xi,yi,Zi,Xr, Vr, Zr, (fik^Ok, tk, Ck) 

Require: Speed of sound c{x) for domain il 

{Algorithm for Shape and Trajectory Reconstruction of 
Moving Obstacles} 

{Estimated time complexity is 0{T^A) where T is the 
number of discretization points for the time of flight, and A 
is the number of discretization points for the angle space} 



for all data points Bk do 

hi ^ ^ 

L = {Xo,Yo, Zq) — {xi,yi,zi) set this initial position to 
be position of transmitter 

S = {aXQ,aYo,aZo) = {xr,yr,Zr) set this initial 
position to be position of receiver 

<^o = <l>k 

©0 = 9k 

To=0 

aTo = 

for s = ^ - 1 do 

{Compute the next point on the ray from the trans- 
mitter by fourth order Runge-Kutta step and the ray 
tracing system 3} 

Xs+i^ RKAxihk,Ts,Xs,Ys,Zs,<i>s,es) 
Ys+i^RKAYihk,Ts,Xs,Ys,Zs,'S>s,es) 
Zs+i = RKAz{hk, Ts,Xs,Ys,Zs,<^s, Qs) 
^s+i=RKU{hk,Ts,Xs,Ys,Zs,^s,Qs) 
es+i=RK4e{hk,Ts,Xs,Ys,Zs,<S>s,es) 
Ts+i ^Ts + hk 

Ps+i = {Xs+i,Ys+i, Zg+i) point on solution of ray 
tracing equations with initial values for transmitter that 
is at time Ts+i away from the transmitter L 
if l(Ps+i G no) then 

There must be a measurment error Continue with 

next data point Bk 
end if 

for all initial angles a$o,a0o in discretized angle 

space of the receiver do 
for p = ^ iV^ - 1 do 

{Compute the next point on the ray from the 
receiver by fourth order Runge-Kutta step and the 
ray tracing system 3} 



aXp-f-i ~ 


- RK4:x{hk, aTp, aXp, aYp, aZp 


,a^p,aQp) 


(lYp+i — 


RK4:Y{hk,aTp, aXp, aYp, aZp, 


a^p, aQp) 


a^p+i = 


RK4:z{hk,aTp, aXp, aYp, aZp, 


a^p,aQp) 


a$p+i = 


RK4:^{hk, aTp, aXp, oYp, aZp, 


a^p, aQp) 


aQp+i ~ 


RK4:e{hk, aTp, aXp, aYp, aZp, 


a<l>p,a9p) 


aTp+i = 


aTp + hk 




P — 


{aXp^i, aYp^i, aZp^i) point on solu- 



tion of ray tracing equations with initial angles 

a$o and a8o and initial position S, that is time 

aTp+i away from S 

if !(P„^^, e flo) then 

Exit this for loop and continue with next pair 
of initial angles a^o, aQo from outer for loop 

end if 

if distance{Ps+i, Pap+i) < ei and jT^+i + 

aTp_|_i ~ tk\ < £2 then 

Pk = Pg+i {Solution for current data point Bk 
found. Continue with next data point Bk+i} 



end if 

if r,+i + aTp+i > tfc + £2 then 

{We are over the travel time budget tk- Con- 
tinue looking for a solution with the next set of 
initial angles a<i>o,a0o.} 
end if 
end for 
end for 
end for 
end for 

RK4 stands for a fourth order Runge-Kutta method although 
other time-dependent numerical methods can be used as well. 
The above algorithm leads to a class of new algorithms 
when data structures are used to store the rays for faster 
processing. For example, when rays from the receiver in its set 
of directions are stored and looked up in a data structure such 
as an array, the performance of the algorithm can be improved 
to 0{TA) where T is the number of discretization points for 
the time of flight, and A is the number of discretization points 
for the angle space. These algorithms are very suitable for 
parallelization and their parallelized implementation is key for 
efficient real-time processing, performance and computational 
efficiency. For example, the computation for tracing the rays 
from the receiver can be parallelized in multiple threads. When 
the above algorithm is run on a set of points {B^} from 
one sampling time interval Tk, the algorithm reconstructs the 
shape of the obstacle during this sampling interval. In order 
to reconstruct the trajectory of the obstacle, the algorithm is 
run on the data points for each of the sampling intervals. The 
method can achieve high resolution because it can process 
a very large number of distinct points on the obstacle's 
boundary. In contrast to tomography where the focus of the 
reconstruction method is to recover the velocity structure of 
the domain, the shape and trajectory reconstruction procedure 
directly finds the shape and trajectory of the obstacle. 

A. Numerical Example 

I have developed a multithreaded Java program that im- 
plements the parallelized algorithm from II using Euler's 
numerical integration method instead of RK4 and reconstructs 
the coordinates of reflection points from input data {Bk} and 
a speed of sound function c{x). Consider a circular reflecting 
obstacle in the plane xy moving in the first quadrant away from 
the origin along the line x=y in a medium with variable speed 
of sound c{x, y) = x + y + 1. We place a transmitter and a re- 
ceiver at the origin and set xl = yl = zl = xr = yr = zr = 0. 
In this case, the domain JIq is a circle of sufficiently large 
radius that contains the origin. Table I shows the computed 
trajectory of a reflection point on the obstacle on the ray path 
from the origin x=y, 1/ > 0, corresponding to data for different 
travel times. 

We check the accuracy of the computation in the above table 
by the Java program as follows. The time for the ray to reach 
to obstacle can be computed by the formula 
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TABLE 1 

Reconstruction of a point moving with speed c{x, y) = x + y + 1 

ON THE line X=Y FOR A FIXED SIGNAL FREQUENCY ^. THE INITIAL 
TRANSMISSION ANGLES ARE </> = f AND 6* = 



Jo c{s) 7o x + y+1 J„ 2x + l 

Therefore, 



By symmetry, for this particular example, this time t is half 
of the total travel time T. Then for a travel time T = 2, or 
t = 1, we compute X = Y = 1.55. This result matches 
the corresponding result for px and py from Table I obtained 
by numerical integration. In the current implementation, the 
reconstructed results are more accurate for shorter travel times 
as can be seen by comparing the result from the table with 
the result from the above formula 8 for T=8 but this accuracy 
can be improved by using numerical integration methods with 
a smaller error 

B. Communication path from receiver to transmitter 

The reconstruction method II above finds the reflection point 
P for the signal from transmitter L to receiver S by finding 
angles of transmission (p and for the receiver at which a 
symmetric signal can be sent to L from S that will reach L 
over the reverse path of the signal detected at S. Therefore, 
the receiver S can communicate with L by sending a response 
signal back to L along angles of transmission (f) and 9. L 
will receive this response if the positions of L or S have not 
changed signigicantly or if the speed of sound has not changed 
significantly in the region containing the communication paths. 
In order to create a more robust connection S can send to L 
several parallel rays that are very close to each other S and L 
can send and receive acknowledgments and create or use an 
existing networking protocol stack. This approach integrates 
easily with current network protocols and architectures. S and 
L recompute the reflection point and transmission angles with 
the receipt of each message. When S and L receive from 
each other control messages with high frequency, S and L 
recompute a current set of transmission angles for routing back 
messages. The frequency of these control messages depends 
on the velocity of S, L and the obstacle, and in the case of 
time-varying speed of sound, on the time derivative of the 
speed of sound or on how quickly the speed of sound changes 
as a function of time. 



III. Network Architecture for High-Performance 
Communications 

Shape and trajectory tracking of moving obstacles that 
reflect ray signals sent from moving Internet hosts enables 
efficient one-hop routing through reflection to the nearest 
Internet router or peer A router in this architecture is a 
transmitter and receiver of ray signals. 

A. Coordinate system 

The architecture is based on a known Cartesian global 
coordinate system. The physical position, for example, the 
high precision GPS coordinates of the transmitter or receiver 
correspond to global coordinate system coordinates xi,yi, zi 
or XrjUr, Zr from the reconstruction algorithms II. 

B. Addressing 

We define a new physical layer networking protocol called 
IHOP. The goal of this protocol is to route messages through 
the environment in at most one hop. The IHOP address of 
a host is always relative to the position of another host, a 
transmitter, and consists of a pair of incident and azimuth 
angles and 6 for reaching the receiver by sending a ray signal 
from transmitter to receiver. The IHOP protocol is defined 
as addressing an IHOP host from an IHOP transmitter by 
the host's IHOP address or by sending messages from the 
transmitter to the host along the direction determined by the 
IHOP angles. 

C. Optimal communication path 

A transmitter communicates with a router via the optimal 
path, the path that takes the least time, for the given speed 
of sound and reflection point. This network architecture is 
illustrated in Fig. I. 

The combination of the methods for obstacle tracking from 
this work with the numerical methods for tomography in the 
presence of obstacles from [L], leads to a new method for 
networking with reflected rays based on tracking the location 
of reflectors and finding the speed of sound in the environment. 
Through reflection, messages can be communicated around 
obstacles blocking the line of sight as well as around areas with 
a high error rate. The reflectors for the shape and trajectory 
tracking method do not have to be installed in order to enable 
a network and can be obstacles such as buildings or cars. The 
two initial angles provide a way of addressing the obstacle and 
the receiver i.e. a way of calling the obstacle and the receiver 
by dialing two angles of transmission. 

IV. Advantages 

A. Improved Security 

The communication path between transmitter and receiver 
is less prone to eavesdropping because communication is via 
a narrow ray. 



B. Increased Network Bandwidth 

The network architecture described in the previous section 
enables high bandwidth networking because a receiver can 
receive information over several rays arriving from different 
directions. Each ray is an independent communications chan- 
nel and delivers information in parallel with other rays detected 
by the receiver. 

C. Increased Reliability, Survivability and Performance 

In order to increase reliability, the transmitter can send 
a message to the receiver along several different rays to 
several different reflection points. The receiver could in this 
case receive the same message from several different rays or 
directions. This provides a natural way for error correction and 
leads to high-performance networking because information 
from messages that are lost or delayed could be delivered 
quickly from a different direction. The fastest message is pro- 
cessed first and this policy improves end-to-end performance. 

D. Reduced delay and better support for high-performance 
networking 

The communication path is computed to match the speed of 
sound in the environment and therefore there will be less en- 
vironmental errors. For example, the effects of temperature or 
humidity differences or other weather conditions are included 
in the computation of the communication path via the function 
for the speed of sound in the environment. The reduced error 
rate will lead to better quality of service. The architecture 
is particularly useful for streaming real-time content to home 
networks or peer-to-peer networks. 

E. Increased Ease of Use 

The new network architecture makes it easier to connect 
home network devices via reflected rays compared to the 
current methods for wireless networking because there is no 
need for additional routers. 

V. Conclusion 

Ray signals may reach a host by an unbroken ray or by a 
broken ray when it reflects from an obstacle. For improved 
bandwidth and reliability, we can reach the host over different 
rays. The IHOP address of the host specifies a ray. 
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Fig. 1. Networking in tlie presence of obstacles. Mobile Internet Phone 1 
transmits a message to Mobile Internet Phone 2 via the broken ray Ray 1. 
Knowing the speed of sound in the environment, Mobile Internet Phone 2 
applies the algorithm for shape and trajectory tracking of moving obstacles 
from II to reconstruct the reflection point at Obstacle 1 and send back response 
signals in the direction of this reflection point. The algorithm from II enables 
the reconstruction of the reflection point as well as the arrival angles of the 
signal at Mobile Internet Phone 2. 
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